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ABSTRACT 

Aims. We present a Gaussianity analysis of the Archeops Cosmic Microwave Background (CMB) temperature anisotropy 
data maps at high resolution to constrain the non-linear coupling parameter /„; characterising well motivated non- 
Gaussian CMB models. We used the data collected by the most sensitive Archeops bolometer at 143 GHz. The instru- 
mental noise was carefully characterised for this bolometer, and for another Archeops bolometer at 143 GHz used for 
comparison. Angular scales from 27 arcmin to 1.8 degrees and a large fraction of the sky, 21%, covering both hemi- 
spheres (avoiding pixels with Galactic latitude |6| < 15 degrees) were considered. 

Methods. The three Minkowski functionals on the sphere evaluated at different thresholds were used to construct a 
X 2 statistics for both the Gaussian and the non-Gaussian CMB models. The Archeops maps were produced with the 
Mirage optimal map- making procedure from processed time ordered data. The analysis is based on simulations of signal 
(Gaussian and non-Gaussian /„; CMB models) and noise which were processed in the time domain using the Archeops 
pipeline and projected on the sky using the Mirage optimal map-making procedure. 

Results. The Archeops data were found to be compatible with Gaussianity after removal of highly noisy pixels at high 
resolution. The non-linear coupling parameter was constrained to f nl = 70+%% at 68% CL and f nl = 70±g!g 5 at 95% 
CL, for realistic non-Gaussian CMB simulations. 
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1. Introduction 

The Cosmic Microwave Background (CMB) radiation is a 
valuable tool to study the early universe. When we ob- 
serve this weak radiation we are observing the universe 
when it was about 400000 years old. Several theories have 
been proposed to explain the origin and the evolution of 
the universe. The Big Bang is the most favored theory 
and it predicts the existence of the CMB. The inflationary 
models are complementary to the Big Bang model. Among 
their many properties we are interested in those related to 
the CMB and the predictions that they prese nt about it . 
The standard, single field, s l ow roll infla t ion ( Guthl Il98lt 
lAlbrecht fe Steinhardtj 119821 : iLindd 1 19821 [[983) is one of 
the most accepted models because of the accuracy of its 
predictions a nd the observations performed with modern 
experiments ()Spergel et al.l [20071 ) . In particular, the stan- 
dard inflation predicts that the primordial density fluctua- 
tions are distributed following a nearly Gaussian distribu- 
tion. These fluctuations are imprinted in the anisotropics 
of the CMB. Therefore if the prediction is correct the 



anisotropies of the CMB should be distributed in a nearly 
Gaussian way. Any detection of non-Gaussian deviations 
would have far-reach ing consequences for our understand- 
ing o f the Universe (|Cruz et al.l l2007ri lYadav fe Wandeltl 
120071 ) . In addition the search for non-Gaussian fluctuations 
in the data is a useful tool to look for residual foreground, 
secondary anisotropies or unidentified systematic errors. 
Some of these effects may introduce non-Gaussian features 
at different levels. 

Many statistical tools have been applied to test the 
Gaussianity of CMB data sets. The Minkowski function- 
als have been applied to different recent experiments 



als have been applied to dinere nt re cent exp eriments 

(ISpereel et all 120071 : Ide Troia et ail 120071 : iGott et alJ[2007t 
ICurto et alJ 120071 ); other examples are the smooth tests 
of goodness-of-fit (ICavon et al.lT2003bt lAliaga et all 120051 : 
iRubino-Martm et al l 120061; ICurto et al l 120071), wavelets 
dBarreiro et al.ll2000t ICayon et alJl2q0ll: IVielva et al.l l2004: 
Mu kheriee fc Wanell2004 ICruz et al.ll2005ll2007al). lo cales- 
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timat ors of the n-point correlations (Eriksen et all 2004. 
1 2005]) . steerable filters dWiaux et all l2006t IVielva et all 
120071 ). and the CMB 1-pdf (jMonteserin et al.ll2007f) among 
others. 
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In this work we analyse the CMB data collected by 
the Archeops experiment. Important results have been ob- 
tained from this experiment since its launc h in 2002. It 
gave t he first link in the Ci determination (jBenoit et al.l 
2003 a[) between t he COBE large angular scale data 
( Smoot et al.lll992D and the first a coustic peak as measured 
by BOOMERanG a nd MAXIMA (|de Bernardis et al.ll2000t 
lHananv et al.|[2000f) . From this it gave a pre cise determi- 
nation of the main cosmological parameters (jBenoit et alj 
2003b])- It also provided an independent confirmation at 
different frequencies of the power spectrum for the range 
£ = 10 to I = 70 of the WMAP first year results 
(jTristram et al.ll2005t ). 

In this study we use the Minkowski functionals 
(iMinkowskil 119031: iGott et al.l 119901 : ISchmalzing fe Gorskil 
11998ft . This new Gaussianity analysis of the A rcheops data 
compl ements the first analysis presented in ICurto et al.l 
( 2007) where only low resolution maps (about 1.8 degrees 
of resolution) were analysed and the /„; constraints were 
imposed on non-Gaussian CMB maps simulated using the 
quadratic Sachs- Wolfe approximation. In contrast, in the 
present work, the constraints on the non-linear coupling 
parameter f n i are obtained using higher resolution (27 ar- 
cmin) Archeops data and realistic non-Gaussian simula- 
tions of the CMB fluctuations wit h the algorithms devel- 
oped bv lLiguori et~aTl (|2003l . l2007l ). 

Our article is organized as follows. Section[2]presents the 
Archeops data and the Gaussian and non-Gaussian simu- 
lations performed. Section 3 describes the statistical meth- 
ods to test Gaussianity and to constrain the f n i parameter. 
In Section 4 we perform an analysis of the instrumental 
noise which will provide the correct masks for our analysis. 
Section 5 is devoted to the /„; constraints and the compar- 
ison with the Sachs- Wolfe approximation. We summarize 
and draw our conclusions in Section 6. 



2. Data and simulations 

2.1. The data 

Archeopfl is a balloon-borne experiment to map the CMB 
anisotropies at high resolution with a large sky cover- 
age. T he experiment is descr i bed i n iBenoit et alj |2002, 
l2004alTbh : IMacias-Perez et al.l (|2007l ). It is based on the 
Planck high frequency instrument (HFI) technology. The 
instrument consists of a 1.5 m telescope pointing at 49 de- 
grees from its vertical axis on board a stratospheric gondola. 
It has 21 bolometers cooled to 100 mK by an 3 He/ 4 He di- 
lution cryostat designed to work at similar conditions to 
the ones expected for Planck. These bolometers operate at 
frequencies of 143, 217, 353 and 545 GHz. After one test 
flight in Trapani (Italy), the instrument flew three times 
from the CNES/Swedish facility of Esrange, near Kiruna 
(Sweden) to Siberia (Russia). The last flight on February 
7 th 2002 provided 12.5h of CMB-quality data for a total 
of 19h. The experiment was performed during the Arctic 
winter to avoid Sun contamination. These data correspond 
to a sky coverage of approximately 30% of the sky, includ- 
ing the Galactic plane. From its four frequency bands the 
two lowest (143 and 217 GHz) were dedicated to the ob- 
servation of the CMB and the others (353 and 545 GHz) 
to the monitoring and calibration of both atmospheric and 
Galactic emission. 

1 http://www.archeops.org 



We used the data collected with two bolometers at 143 
GHz. We used the data of the 143X03 bolometer for the 
main analysis (see Sect. [5]) and the data of the 143X03 
and 143X04 bolometers for the noise analysis (see Sect.|4|). 
However, for the /„; constraints, the data from the sec- 
ond bolometer 143X04 are not used due to their higher 
noise level and worse systematic erro rs. The characteris- 
tics of the bolometers are des cribed in IMacias-Perez et ail 
(|2007l ). In ICurto et all (j2007f ) the Gaussianity analysis of 
Archeops data was performed at lo w resolution, in partic- 
ular HEALPix (jGorski et alj I2005D N stde = 32. Here we 
complement that work and analyse the data at the same 
and higher resolutions: N S id e — 32, 64, and 128 using the 
realistic non-Gaussian simulations presented below. 

First we processed and cleaned the Time Ordered 
Inform ation (hereafter TOI) as described in lTristram et al.l 
(2005?). Then the data maps at different resolutions were 
produced using the M irage optimal map-making procedure 
(|Yvon fc~M avct 200jf). All the analyses presented here were 
performed on a fraction of the Archeops observed region af- 
ter masking out pixels with Galactic latitude |6| < 15°. This 
corresponds to 21% of the total sky. Unlike the analysis in 
ICurto et al.l (|2007| ). restricted to north Galactic latitudes, 
south Galactic latitudes are included in the analysis. 

2.2. Gaussian simulations 

We have performed Gaussian simulations of t he CMB sig- 
nal an d of the Archeops noise as described in ICurto et alj 
(2007) . The noise simulations were obtained from Gaussian 
random realisations of the time-domain noise power spec- 
trum of the bolometers. The constructed noise time ordered 
data were then projected using the Mirage optimal map- 
making procedure. The CMB Gaussian simulations were 
obtained from random realisations of the best-fit Archeops 
CMB power spectrum. From these maps we constructed 
Archeops time ordered data and then we project them us- 
ing Mirage. In this way, the filtering of the Archeops data 
was taken into account. 



2.3. Non-Gaussian CMB simulations 

There are several alternative theories to the standard in- 
flation theory that introduce non-Gaussian fluctuations in 
the CMB. One simple model that describes weakly non- 
Gaussian fluctuations in matter and radiation is obtained 
by introducing a quadratic term in the primordial gravita- 
tional potential (ISalopek fc Bondlll990l:lGangui et al.|[l994j: 
IVerde et al.l 120001 : iKomatsu fc SpergelfeoOlh 



*(x) = $ L (x) + f n iM(x) - <$!(*)}} 



(1) 



where <&l{x) is a linear random field which is Gaussian dis- 
tributed and has zero mean, and f n i is is the non-linear cou- 
pling parameter. Non-Gaussianity of this type is generated 
in v arious classes of non -standard inflationary models (see, 
e. g. iBartolo et aT]|2004 for a review). To obtain the CMB 
anisotropies generated by such a primordial gravitational 
poten tial we use the algorithm descr ibed in ILiguori et al.l 
(|2003l ) for the temperature maps or ILiguori et al.l ( 20071 ) 
for temperature and polarization maps. 

In this case the multipole coefficients ai m of the CMB 
temperature map can be written as 
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where ajf? is the Gaussian contribution and a[^[ G ' is the 
non-Gaussian contribution. 

For our simulations we use the ACDM model that 
best fits the WM AP data and a modif i ed ve rsion of the 
CMBFAST code (jSeliak fc Zaldarriagal Il996fl to obtain 
the Gaussian and non-Gaussian contributions as described 
above. We produced a set of 300 high resolution full sky 
temperature Gaussian maps ATq and their complementary 
non-Gaussian maps AT/vg a t the same resolution. The total 
CMB map with a non-Gaussian contribution is therefore 

AT = AT G + f ni AT NG . (3) 

To constrain /„; with Archeops data we need to transform 
these simulations into Archeops simulations. For this pur- 
pose, each simulated non-Gaussian map given by Eq. [3] 
was converted into Archeops time-ordered data account- 
ing for the Archeops pointing. These time-ordered data 
were re-normalised to account for intercalibration errors 
between Archeops and WMAF0, and added to Gaussian 
simulations of the Archeops instrumental noise constructed 
as above. Then the time ordered data was projected on 
the sky using Mirage. Although, the Mirage algorithm is 
weakly non- linear, i.e., the map obtained from ATg cuid 
ATng separately is slightly different from the one obtained 
from the sum of the two contributions, we have tested that 
the Minkowski functionals are very similar in both cases. 
Thus, for our analysis, we can construct the Archeops sim- 
ulations for different values of /„; by adding the Mirage 
map for ATq to the one for AT^g multiplied by /„;. This 
corresponds to an important saving of CPU time as the 
Mirage algorithm is a significantly time consuming process 
(|Curto et al.ll2007fl . 

Thus, we transformed our sets of 300 ATg and AT/vg 
simulations into CMB Archeops simulations (s g and s ng ) 
at the three considered resolutions N S ide =32, 64, and 128. 
The final non-Gaussian simulations accounting for f n i were 
computed from s g + f n is n g + % where n g corresponds to 
the Archeops Gaussian noise simulations. 

3. Methodology 

The Gaussianity analysis and the constraints on the f n i 
parameter were performed using the Minkowski function- 
als. Detailed theoretical information about thes e quan - 
tities is presented for ex ample in iGott et aTl |l990); 
ISchmalzing & Gorskil (l998). Some recent examples of ap- 
plications of the M inkowski functionals to the CMB are in 
ICurto et alj (l2007f) fo r Arch eops data at low angular reso- 
luti on. Ide Troia et all (120071) for BOOMERanG 2003 data, 
and lSpergel et alj (|2007l )for WMAP 3rd year data. 

For a scalar field AT(n) in the sphere we have three 
Minkowski functionals given a threshold v. Considering the 
excursion set of points Q v where AT in) j a > v the three 
Minkowski functionals are: the area A(v) of the contour 
length C(v) of Q V: and the genus G(v) (proportional to the 
difference between hot spots above v and cold spots below 
v). The expected values of the functionals fo r a Gaussian 
random field are (ISchmalzing fc G6rskilfl998h 

(My)) = - f 1 --= y o exp(-t 2 )dt) 

2 The normalisation is simply obtained by multiplying the 
si mulated TOI by a co nstant factor / = 1/1.07 as described 
in iTristram et~alTl|2005l ). 
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where a = J2}ZT ( 2e + l ) C t and r = T/tZT ( 2e + + 
l)/2. These expressions are only valid in an ideal case. 
In practice, we will obtain these quantities from Gaussian 
CMB simulations in order to take into account the noise, 
mask and pixel effects. 

For a Gaussian random field the Minkowski functionals 
are approximately Gaussian distributed, therefore we can 
use a x 2 statistic to test the Gaussianity of a CMB map. 
This statistic is computed with the three Minkowski func- 
tionals evaluated at n t h different thresholds. 

3n t h 

x 2 = £K-(^>)^>i-<^->) (5) 

where () is the expected value for the Gaussian case, v is a 
3n t ^ vector 

(Afa), ..,A(v nth ), C{v x ), .., C{v nt J, G{v x ), .., G{v nth )), (6) 
and C is a 3n t h covariance matrix 

Cij = (viVj) - (Vi){vj). (7) 

The Gaussianity analysis consists of computing the \ 2 
statistic of the data and compare it with the value of the \ 2 
statistic of a set of Gaussian simulations of the data (CMB 
signal plus instrumental noise). 

Another important analysis is the estimation of the /„; 
parameter. In this case, we can use a \ 2 test with the 
Minkowski functionals 

3n th 

X 2 (fnl) = J2 ^ " (^fn^tfnlXVj ~ (Vj)f nl ) (8) 

i.i=l 

where () f nl is the expected value for a model with /„;, and 



Cij {fni 



> J * V 3lU 



)f nl . For low values of f n i 



(fni < 1500) we have that Cij{f n i) ~ Cy(/„i = 0) = C!y, 
and therefore we can use the approximation 

X 2 (fnl) = £>i- " (»*>/»,)• (9) 

The best-fit f n i for the Archeops data is obtained by mini- 
mization of X 2 (fni)- Error bars for it at diffirent confidence 
levels are computed using the Gaussian simulations. 

Other than the analysis of the data at each resolution 
separately we can also analyse the data by combining the 
information at different resolutions. Assuming maps at n res 
different resolutions we can define a vector 

V = (vi,V 2 , —,Vn re .), (10) 

where each Vi corresponds to the vector given by Eq. [5] for 
the resolution i. With this combined vector we can compute 
X 2 statistics in the same way as in Eqs.[S]and|n] 



3N th 

X 2 = YAVi-mCrMVj-M)) 



X 2 (fnl 



3N th 

= Ew 



(11) 



(12) 
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where () is the expected value for the Gaussian case, () f nl is 
the expected value for a model with f n i , Nth — J2k=i n th ' 

(k) 

n th is the number of thresholds used at resolution k, and 

Cij = (v i v j )-W(v j ). 

4. Gaussianity analysis of the Archeops 
instrumental noise 

From M arias-Perez et al.l (|200?T) it is expected that the 
noise at a given pixel behaves as Gaussian when the pixel 
has been observed a sufficiently large number of times, 
about 100 times. In order to check the noise behaviour, 
we perform a Gaussianity analysis of the Archeops instru- 
mental noise. This is needed to exclude from the main 
Gaussianity analysis the possible scales and regions of the 
sky where the Archeops noise is non-Gaussian. 

4.1. Two bolometer analysis 

A simple way to analyse the instrumental noise is to sub- 
tract the data of two different bolometers at the same 
frequency so that the sky signal is removed in the final 
map. We used the data collected by two different Archeops 
bolometers at 143 GHz: 143K03 and 143K04. As the CMB 
contribution is the same for both bolometers but the noise 
is not, after subtraction and assuming no systematic errors 
only the noise contributions remain. Defining dxo3 as the 
data collected by 143K03 and dxo4 the data collected by 
143K04, the map that we analysed is 

d = dx 03 — dxoA ~ n-K03 — n>K04- (13) 

We analysed the difference data maps given by Eq. [T3] for 
three different resolutions: N S id e = 32, N S id e = 64, and 
Nside = 128. To construct simulations of the difference 
map, we used a set of 10 3 Gaussian signal simulations and 
10 3 Gaussian noise simulations for each bolometer and each 
resolution. These simulations are used to compute the mean 
value (vi) and the covariance matrix Cy of the Minkowki 
functionals as described in Sect. [3] 

In Fig. Q] we present, from left to right, the Gaussianity 
analysis of the difference maps at the three resolution con- 
sidered, N side = 32, N slde = 64, and = 128. The 
histograms correspond to the values of the % 2 obtained for 
10 3 simulations of difference maps. The vertical lines corre- 
spond to the value of \ 2 for the Archeops data. The solid 
lines are the theoretical vi„ distribution normalised to 
the number of simulations and the size of the binned cell. 

From this analysis we can see that the data given by 
Eq. [13] becomes non-Gaussian at high resolution. This is 
most probably due to highly noisy pixels in the difference 
maps corresponding to regions of the sky observed with 
little redundancy. A previous analysis (jMacias-Perez et alj 
l2007f ) has shown that the Archeops instrumental noise in 
the map domain is Gaussian distributed for a given pixel 
when this pixel has been observed a significant amount of 
time, typically above a few hundred independent observa- 
tions (hits) per pixel, although no precise estimate of the 
required number of hits per pixel was given. 

This can be easily done using the statistical tools pre- 
sented in this paper. For this purpose we perform our anal- 
ysis excluding highly noisy pixels defined as those pixels 
presenting a number of hits below a given threshold. We 



computed the % 2 statistic of the data (Eq. [TB"]) , its cumu- 
lative probability, and the remaining area using different 
thresholds for the number of hits. By comparing to the 
Gaussian simulations we observe that for both N S ide = 64 
and N S id e = 128, increasing the threshold of the num- 
ber of hits implies that the data becomes compatible with 
Gaussianity. In particular, for N S ide = 64 the data become 
compatible with Gaussianity when pixels with fewer than 
250 hits are removed (leaving 46% of the original area) and 
for N S ide = 128 the data start to become compatible with 
Gaussianity when pixels with fewer than 150 hits are re- 
moved (leaving 13% of the original area). 

4.2. Single bolometer analysis 

From the above analysis we have proved that highly 
noisy pixels are responsible for the non-Gaussianity of the 
Archeops noise at high resolution. In order to obtain proper 
limits for the /„; parameter we need to exclude these pixels 
in our final analysis of the 143X03 map. 

To avoid any contamination from the 143X04 noise map 
when characterizing the noise map for the 143K03 bolome- 
ter the latter was computed using the WMAP data as fol- 
lows 

n K 03 = d K Q 3 - WM AP K03 * / (14) 

where WMAP kqs is the combined WMAP data computed 
as described in ICurto et al.l (|2007t ) and converted into an 
Archeops map as described in Sect. 12.31 / is an intercalibra- 
tion param eter between WMAP a nd Archeops, / = 1/1.07, 
as given bv lTristram et alj (|2005h . In this the noise 

on the WMAP data at the considered map resolutions is 
negligible compared to the Archeops noise, we expect no 
contamination. 

We analysed the noise map given by Eq. [TJ] for three 
different pixel resolutions: N S ide — 32, N S id e — 64, and 
N S ide = 128. For each resolution, we also analysed a set of 
10 3 Archeops noise simulations uko3- The results are pre- 
sented in Fig. [2] The histograms correspond to the \ 2 °f the 
Minkowski functionals of a set of 10 3 nxo3 Gaussian simu- 
lations, vertical lines correspond to the data maps given by 
Eq. [14] and the solid lines are the Xanth distributions. We 
can see that for low resolution (N S ide = 32 and N S ide = 64) 
the noise map is compatible with the Gaussian noise simu- 
lations. Therefore for these two resolutions we can use the 
full available area in the f n i analysis. 

In the case of N S ide — 128 the noise map is not compat- 
ible with Gaussianity. As above, we reanalysed this noise 
map removing highly noisy pixels with a number of hits be- 
low a given threshold. This analsysis is presented in Fig. [3] 
where from left to right we plot the % 2 statistic of the 
143K03 noise map, its cumulative probability, and the frac- 
tion of available area for different values of the minimum 
number of hits. We can see that the Archeops noise for 
the 143AT03 bolometer becomes clearly Gaussian when we 
remove pixels with a number of hits lower than 90. This 
means that the area where the noise is Gaussian is 32% of 
the original area at N S ide = 128. In conclusion, the regions 
of the sky that we can use for the Gaussianity analysis are 
the whole initial area at N s ide = 32 and N S ide = 64 and 
32% of the initial area at N S ide = 128. 

The 143AT04 bolometer has a qualitatively similar be- 
haviour for the noise analysis. That is, the noise map is 
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Fig. 1. Distribution of the x 2 values from the Minkowski Gaussianity analysis for the Archeops 143K03-143K04 map at 
N S ide = 32, N S ide = 64, and N S id e = 128. The vertical lines show data, the histograms are obtained from a set of 10 3 
simulations, and the solid lines are the expected x 2 distributions for 3rith degrees of freedom. 
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Fig. 2. Distribution of the x 2 values from the Minkowski Gaussianity analysis for the Archeops 143K03 noise map at 
N S ide = 32, N si de — 64, and N si d e — 128. Vertical lines show data, the histograms corresponds to sets of 10 3 noise 
simulations, and the solid lines are the x 2 distribution with 3n t h degrees of freedom. 




Fig. 3. From left to right, \ 2 statistic of the 143K03 noise map, its cumulative probability, and the fraction of available 
area for different values of the minimum number of hits at N S id e = 128. This map becomes compatible with Gaussianity 
when the minimum number of hits is approximately 90. 



compatible with Gaussianity at low resolution and is non- 
Gaussian at high resolution. Nevertheless as the noise for 
this bolometer has more systematic errors than the 143if03 
bolometer, the thresholds for the number of hits where 
the maps become compatible with Gaussianity analysis are 
higher. This implies that the available area where the noise 
is Gaussian is smaller. 



5. Results 

In this section, we perform the Gaussianity analysis of the 
143X03 Archeops bolometer data in order to constrain 
the non-linear coupling parameter, /„;. For this, we first 
consider realistic non-Gaussian simulations as described in 
Sect. 12.31 Then the results are compared to the ones ob- 
tained for non-Gaussian simulations for whi ch only the 
Sachs- Wolfe contribution is included, as in ICurto et alj 
(|2007f ). 



5.1. Gaussianity analysis 

For this analysis we apply the statistical methods described 
in Sect. [3] to the Archeops 143K03 bolometer data. Figured] 
shows the Minkowski functionals as a function of threshold 
for maps with N side = 32, N side = 64, and N side = 128. 
In addition, we plot the acceptance intervals for the 68% 
(inner), the 95% (middle), and 99% (outer) significance lev- 
els obtained from 10000 noise and signal Gaussian simula- 
tions of the Archeops 143K03 data as described in Sect. O 
For all the analysed cases the data are compatible with 
Gaussianity at least at the 99% significance level. 

In Fig. Owe present the results of the Gaussianity anal- 
ysis of the Archeops 143K03 bolometer data at different 
resolutions using the % 2 test described by Eqs.lBlandfTTl In 
particular, we analysed the data at N S id e — 32 (11 thresh- 
olds, 3 functionals), N si d e = 64 (13 thresholds, 3 function- 
als), N s ide = 128 (15 thresholds, 3 functionals), the com- 
bination of the data at N S id e = 32 and N S id e = 64 (to- 
tal of 24 thresholds, 3 functionals), and the combination 
of the data at N side = 32, N slde = 64, and N stde = 128 
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Fig. 4. From left to right, and from top to bottom, the area, contour length, and genus of the data (asterisk *) as a 
function of threshold for the 143K03 maps with N S id e = 32, N S ide = 64, and N s id e = 128. We also plot the acceptance 
intervals for the 68% (inner in red), the 95% (middle in green), and 99% (outer in magenta) significance levels given by 
10000 Gaussian simulations of signal and noise. 
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Fig. 5. From left to right, and from top to bottom, the Gaussianity analysis of Archeops data maps at resolutions: 
Nside = 32, N S id e = 64, N s ide = 128, the combinations N S id e — 32,64, and the combinations = 32,64, 128. The 

histograms correspond to the x 2 of 2000 Gaussian simulations, the vertical lines are the \ 2 of the data and the solid lines 
are the expected \ 2 distribution with 3N t h degrees of freedom. 



(total of 39 thresholds, 3 functionals). To avoid confusion 
by the non-Gaussianity of the Archeops noise at high res- 
olutions, the high noise pixels were excluded as discussed 
in the previous Section. The histograms in the plot corre- 
spond to the analysis of 2000 Archeops signal plus noise 



Gaussian simulations at each resolution, the solid lines cor- 
respond to the expected Xs\n distribution, and the vertical 
lines correspond to \ 2 values of the data at each resolution. 
For all the resolutions and combination of resolutions, the 
Archeops data are consistent with Gaussianity as expected 
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Table 1. \ 2 f° r the three Minkowski functionals for different resolutions and thresholds. 



Resolution 


Area (%) 


2 

Xdata 


Degrees Of Freedom 


(x 2 ) 


V((x 2 - (x 2 )) 2 ) 


P(X 2 < Xdata) 


32 


100 


43.62 


33 


33.55 


10.05 


0.868 


64 


100 


44.39 


39 


40.55 


10.33 


0.679 


128 


32 


50.59 


45 


46.29 


10.79 


0.688 


32,64 


100,100 


85.14 


72 


73.88 


14.29 


0.804 


32,64,128 


100,100,32 


134.05 


117 


120.30 


18.40 


0.791 



Table 2. Best fit /„; of Archeops data at different resolutions, mean, dispersion and some percentiles of the /„; distri- 
butions obtained from 2000 Gaussian simulations. 



Resolution 


Area (%) 


Best Fit f nl 


(fnl) 




^"o.ieo 


^0.840 


Ao.025 


Xo.975 


32 


100 


90 


48 


702 


-670 


785 


-1370 


1445 


64 


100 


75 


-29 


645 


-665 


615 


-1315 


1270 


128 


32 


-45 


112 


880 


-810 


1000 


-1665 


1805 


32,64 


100,100 


70 


-3 


550 


-525 


535 


-1130 


1070 


32,64,128 


100,100,32 


70 


19 


503 


-470 


520 


-990 


1005 




Fig. 6. From left to right, X 2 (fni) of Archeops data vs f n i and the best fit value of f n i for a set of 2000 Archeops Gaussian 
simulations (sfco3 + ^fc03) combining the maps at resolutions N si d e — 32, 64, 128. 



from the previous figure. In Table[T]we present the \ 2 value 
of the data for the different cases that we have analysed, 
the total number of degrees of freedom, the mean and the 
dispersion of the \ 2 value of the Gaussian simulations, and 
the cumulative probability of the data computed from the 
distribution given by the simulations. 

If we use all the available area at N S id e = 128 we find 
that Archeops data are not compatible with Gaussianity. 
In particular, we have obtained P(x 2 < Xdata) — 0.995 
analysing the data at N slde = 128, and P(x < Xdata) = 
0.994 analysing the data at N S i de — 32, 64, 128. These non- 
Gaussian deviations can be clearly associated with the non- 
Gaussianity of the Archeops noise at high resolution due to 
highly noisy pixels. 

5.2. Constraints on f n i for realistic non-Gaussian simulations 

The constraints in the f n i parameter were obtained using 
300 CMB non-Gaussian simulations such as the one de- 
scribed in Sect. 12.31 and applying to the Archeops 143K03 
bolometer data the % 2 test in the case of non-Gaussian fluc- 
tuations defined by an /„/ parameter (Eq.[9]). We computed 
an Archeops TOI for each realization and from it we com- 
puted the corresponding Mirage Archeops simulation for 
the 143X03 bolometer. We did this for the Gaussian and 
non-Gaussian part separately. So we have 

dkmifnl) — S g ,fe03 + fnl * S ng ,k03 + "-g,fe03 (15) 



where s g ,K03 is a Gaussian CMB simulation, s ng .K03 is its 
corresponding non-Gaussian part, and n gt K03 is a Gaussian 
instrumental noise simulation. 

We computed the mean value of the Minkowki func- 
tionals, (V)f nl , for -1500 < f n i < 1500. We assumed 
that in this interval the covariance matrix associated with 
them was dominated by the Gaussian contribution: i.e. 
Cij(fni) — Cij(f n i — 0). Therefore, it was computed 
from 10 4 Gaussian simulations of signal and noise (s g} k03 + 
^3,fe03) of the Archeops 143K03 bolometer. We obtained 
the X 2 (fm) of Archeops data given by Eqs.[51and[T2lfor the 
same combination of resolutions as the ones described in the 
above Gaussianity analysis. In all cases we find the value of 
f n i that minimizes % 2 (fnl)- This is the best-fit value for the 
f n i parameter. The significance of these values is estimated 
using 2000 Gaussian simulations. For each simulation we 
computed x 2 (fni) vs f n i and obtained its best fit f n i- At 
the end we have a set of 2000 values of this parameter. 
As the simulations are Gaussian, these values are centred 
around f n i = 0. 

Table [2] lists the best-fit f n i to the data for each case 
analysed. We also present the main properties of the distri- 
bution of the fni parameter as obtained from the simula- 
tions. The smaller dispersion corresponds to the case where 
we combine the data at the three resolutions N S id e =32, 64, 
and 128 (only 32% of the available area). Therefore this case 
leads to the best constraints of f n i for the Archeops data. In 
Fig. [HI we present the X 2 (fni) vs /„; of the Archeops data 
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and the histogram of the best fit /„; obtained from 2000 
Gaussian simulations for this optimal case. From this we 
conclude that f nl = 70t^ at 68% CL and f nt = 70±<% 5 
at 95% CL. 

We may wonder if the constraints on the f n i parameter 
could be improved by increasing the area at N si d e = 128. 
We have shown that, if the noise were Gaussian, including 
the whole area available would have produced very similar 
constraints. The reason is that the excluded pixels were 
the noisiest ones and therefore increasing the area would 
not improve the results. 

5.3. Comparing with the Sachs-Wolfe approximation 

On large angular scales the main contribution to the 
CMB anisotropi cs is given by the Sachs Wolfe effect: 
AT/T = -<j>/3 (|Sachs fc Wolfdll967t iKomatsu fe Spergell 
l200lh . where </> is the primordial potential. Therefore, in this 
case, the temperature map can be approximated in terms 
of a linear and a non-linear contribution in the following 
simple way: 

where AT^ is Gaussian. We call this the Sachs- Wolfe ap- 
proximation and to avoid confusion hereafter the f n i pa- 
rameter in this case is called . Thi s approximation has 
been used in several works (for examplelCavon et alj | 2003at 
ICurto et al.ll2007l: Ueong fc SmootJl2007l) . For this approxi- 
mation the Archeops non-Gaussian simulations can be ob- 
tained as follows 

3 f sw 

403(/«0 = S 3: k03 (4,k03 ~ ( s s,fe03>) + n g,k03, (17) 

where s 3t k03 an d n g ^ko3 are Gaussian CMB and noise sim- 
ulations respectively. As discussed above, this approxima- 
tion is only valid for angu lar resolutions of a few degrees 
(|Komatsu fc Spergell2001[ ) . Nevertheless we performed the 
non-Gaussianity analysis of the Archeops data at all con- 
sidered scales in order to constrain and then tested 
the validity of this approximation as a function of the 
scale. In Table [3] we present the results of this analysis. 
We can see that for all the considered resolutions, the /„/ 
constraints are better than for the case with the realistic 
simulations. This difference becomes more important as the 
resolution increases. The best constraints on /„; in this case 
are = 25±?°° at 68% CL and f^ w = 25±^ at 95% 
CL using the combination of the data at N S id e = 32, 64 , 
and 128 (only 32% of the available area). 

In general the Sachs- Wolfe approximation overestimates 
the non-Gaussianity of the CMB fluctuations for a given 
/„/ value. This can be clearly seen in Fig. [7] where we 
present from left to right the histogram of the non-Gaussian 
part of a full sky CMB realistic simulation ATjvg and the 
corresponding quadratic approximation — 3/„i/To(ATg. — 
(ATq) 2 ) at resolution N si d e = 16 i.e. a pixel size of 3.6 
degrees. We can see that even at this scale the simulations 
are very different. The approximation is just the square of 
a Gaussian distribution centred to have zero mean whereas 
the exact simulation is more Gaussian-like. This explains 
why the constraints on the /„; parameters are tighter for 
the Sachs-Wolfe approximation than for the exact case. 



Therefore, the results obtained with the Sachs- Wolfe ap- 
proximation should be taken with caution, since the error 
bars are clearly underestimated. 

6. Conclusions 

In this paper we have presented a complete non- 
Gaussianity analysis of the Archeops data at 143 GHz us- 
ing the Minkowski functionals. First, we characterised the 
Archeops instrumental noise by taking the difference of the 
data of the two most sensitive bolometers, 143K03 and 
143K04. From this we found non-Gaussian deviations at 
high resolution, 27 arcmin (N S id e = 128 in the Healpix pix- 
elization scheme). This is due to the noisiest pixels for which 
at high resolution the number of observations per pixel 
does not allow good systematic error removal. A more de- 
tailed analysis has been performed for the 143i4T03 bolome- 
ter for which a noise map was obtained by subtracting the 
WMAP CMB data. From this analysis we found that pixels 
with a number of observations below 90 were non-Gaussian. 
Masking out those pixels, the noise map is compatible with 
a Gaussian model. Similar results were obtained with the 
143-KT04 bolometer although the minimum number of ob- 
servations per pixel for Gaussianity was of the order of 150 
and was not considered for further analysis. 

Second, masking out these highly noisy pixels we per- 
formed a Gaussianity analysis of the Archeops 143iT03 data 
at low and high resolution. We found that the data are 
compatible with Gaussianity at N S id e — 32, N S id e = 64, 
N S ide = 128, and for the combinations N S id e = 32, 64, and 
for the combinations N S id e = 32, 64, 128 at better than 95% 
CL. From t his analysis a n d usi ng realistic non-Gaussian 
simulations (|Liguori et al.l (|2003f )) we imposed constraints 
on the /„; parameter at these resolutions. The tightest con- 
straints are f nl = 70l^° at 68% CL and f nl = 70±l° 2 7 5 at 
95% CL. 

Third, we also imposed constraints on f n i using the 
Sachs- Wolfe approximation, = 25^°° at 68% CL and 

= 251^00 at CL. For comparison notice that these 
constraints are a f actor of k 3 smaller than those given in 
ICurto et alj (2007) where only low resolution (N S id e — 32) 
maps were considered. 

Finally, we performed a detailed comparison of the real- 
istic non-Gaussian simulations used in this paper and those 
from the Sachs- Wolfe approximation. From this we con- 
clude that even at low resolution the Sachs- Wolfe approxi- 
mation overestimates the non-Gaussianity of the CMB fluc- 
tuations and therefore the f n i constraints imposed are too 
tight by a factor of three, as shown above. 
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Table 3. Best fit of Sachs- Wolfe /„; of Archeops data at different resolutions, mean, dispersion and some percentiles. 



Resolution 


Area (%) 


Best Fit f ni 


(fnl) 


v / </^> - (Ui) 2 


^0.160 


^0.840 


^0.025 


-^0.975 


32 


100 


-150 


33 


418 


-375 


425 


-800 


875 


64 


100 


400 


-3 


252 


-250 


250 


-500 


500 


128 


32 


-175 


10 


232 


-225 


250 


-450 


450 


32,64 


100,100 


200 


10 
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-225 
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500 


32,64,128 


100,100,32 


25 


8 


170 


-175 


175 


-325 


350 




-2x10 



-1x1 0"° O 
T(K) 



1 x1 CT 



1 500 



1 OOO - 




500 - 



1x10 



1x10" 



T(K) 



Fig. 7. From left to right, the non-Gaussian part of a full sky CMB simulation at low resolution (N S id e — 16) and the 
corresponding Sachs- Wolfe quadratic approximation. 
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